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ABSTRACT 

In this paper we investigate a hitherto unexplored source of potentially significant error in 
stellar population synthesis (SPS) models, caused by systematic uncertainties associated with the 
three fundamental stellar atmospheric parameters; effective temperature Tf.ff, surface gravity g, 
and iron abundance [Fe/H]. All SPS models rely on calibrations of Tf,ff, logg and [Fe/H] scales, 
which are implicit in stellar models, isochrones and synthetic spectra, and are explicitly adopted 
for empirical spectral libraries. We assess the effect of a mismatch in scales between isochrones 
and spectral libraries (the two key components of SPS models) and quantify the effects on 23 
commonly used diagnostic line indices. We find that typical systematic offsets of lOOK in Te//, 
0.15 dex in [Fe/H] and/or 0.25 dex in \ogg significantly alter inferred absolute ages of simple 
stellar populations (SSPs) and that in some circumstances, relative ages also change. Offsets 
in Teff, logg and [Fe/H] scales for a scaled-solar SSP produce deviations from the model which 
can mimic the effects of altering abundance ratios to non-scaled-solar chemical compositions, 
and could also be spuriously interpreted as evidence for a more complex population, especially 
when multiple-index or fuU-SED fitting methods are used. We stress that the behavior we find 
can potentially affect any SPS models, whether using full integrated spectra or fitting functions 
to determine line strengths. We present measured offsets in 23 diagnostic line indices and urge 
caution in the over-interpretation of line-index data for stellar populations. 

Subject headings: stars: evolution — galaxies: evolution — galaxies: stellar content 



1. Introduction 

In recent years, stellar population synthesis 
(SPS) models have become a fundamental tool in 
the study of both Galactic (resolved) and extra- 
galactic (unresolved) stellar populations. Pre- 
dictions of both photometric and spectroscopic 
properties, in terms of broad-band colors and 
line index strengths, from various SPS models 
have yielded methods for tackling the well known 
age-metallicity degeneracy in simple (single-age, 
single n ietallicity) stellar populations (SSPs - 
see e.g. iPercival et al.ll2009t iKotulla et al.|[200g : 



Coelho et all 120071: IVazguez fc Leithereij 12005: 
Marastonll2005t iBruzual fc Charlotll2003l for 
recent examples). 



some 



As observational data, and in particular spec- 
troscopic data, are becoming more precise, the 
models are becoming more sophisticated and 
methods are being developed to determine more 
detailed information such as the level of a- 
enhancement in SSPs, chemical evolution, and the 
general properties of the star formation histories of 
composite populations. However, results from dif- 
ferent sets of models often disagree and there are 
currently many unresolved issues concerni ng their 
i nterp retation and implementation. Carte r et al 



( 20091 ) investigated the ability of several "off the 
shelf" SPS models to reproduce optical and near 
infra-red colors (from u through to K) of a small 
sample of well-studied nearby elliptical galaxies. 
Fitting to SSPs, making the implicit assumption 
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that giant ellipticals can be well represented by 
a simple population with no significant star for- 
mation history, they went on to compare the ages 
and metallicities predicted from the various mod- 
els, using a simple reduced test. Carter et al. 
found that, for a similar x^, best-fits from dif- 
ferent SPS models can give ages which differ by 
5 Gyr or more whilst predicted [Fe/H] can easily 
differ by 0.4 dex (a problem here being the [Fe/H] 
sampling, which for most sets of SPS models is 
limited to ~0.3 dex). Predictions from spectro- 
scopic data, in the form of diagnostic line indices, 
can have similar (if not worse) discrepancies, as 
demonstrated by the results of the lAU Sympo- 
sium 241 Stellar Population Challcng({3- For the 
well-studied Galactic Globular Cluster 47 Tuc, 
which was set as a test o bject with an o bserved 
spectrum provided from ISchiavon et al ] ioos), 
predicted ages from various models range from 
5 Gyr to >18 Gyr whilst predicted [Fe/H] ranges 
from —0.1 to —0.8 dex. 

Several recent studies have focussed on discrep- 
ancies that can arise because of properties of the 
observed population, which may not be fully ac- 
counted for in the models. These include ex- 
tended blue horizontal branches and blue strag- 
glers which ca n both make a population look spu- 



rious ly young (jSchiavon et al.ll2004t iCenarro et al.l 
2008L respectively, using diagnostic line indices), 



and thermally pulsing AGB stars which strongly 
affect the near-infrared flux (and hence colors) 
at intermediate ages, significantly altering the in- 
ferred galaxy masses (Bruzual 2007). However it 
is clear that there can be significant offsets be- 
tween ages and metallicities determined from (or 
predicted by) different sets of SPS models which 
must, at least in part, be due to systematic dif- 
ferences in the model ingredients. This problem 
is clearly demonstrated by the 'zero-point' offset 
in some models utilizing diagnostic line indices 
which means that inferred ages of Galactic glob- 
ular clusters can be older than 14 Gyr. Although 
these pr oblems have been noted by several groups 
(seee.g.LV azdekis et al.| l200ll : ISchiavon et al.l[20o3 



Cenarro e t al. 2008) little work has been done to 



determine the specific sources of these offsets or 
quantify their effects in terms of the true system- 
atic errors which should be assigned to derived 



-'http: / /www. astro. rug.nl/~sctrager/challenge/ 



ages and metallicities. 

Physical parameters of stars within a popula- 
tion, such as effective temperature, Te//, surface 
gravity, logg, and iron abundance, [Fe/H], are key 
parameters in the construction of the SPS models 
(see Section[2]) and one, as yet, unexplored aspect 
of SPS is the adopted temperature and metallic- 
ity scales used in the underlying stellar models and 
spectra. When generating stellar evolution mod- 
els (which give rise to isochrones, a fundamental 
building block in SPS modelling), the temperature 
and metallicity scales effectively come out of the 
models themselves. This is also the case for syn- 
thetically generated stellar spectra, for which stel- 
lar atmosphere models and adopted line lists give 
rise to some intrinsic temperature and metallicity 
scale. For empirical spectra, a specific tempera- 
ture and metallicity scale must be adopted in or- 
der to determine physical parameters of the stars. 
There can potentially be mismatches between the 
temperature and metallicity scales of the under- 
lying stellar models and those of the stellar spec- 
tra (whether synthetic or empirical), but any im- 
pact this may have on the resulting SPS models 
has not yet been investigated systematically and 
quantified. The aim of this paper is to address this 
issue. 

The structure of the paper is as follows: In Sec- 
tion [2] we review the ingredients required to con- 
struct SPS models whilst Section [3] describes some 
preliminary tests which motivated this work. Sec- 
tions m and [5] describe the specific tests applied 
to the models and the results of these tests whilst 
Section [6] briefly assesses the significance of the re- 
sults in the context of typical observational errors 
for extragalactic systems. In Section[7]we summa- 
rize our results and discuss their implications for 
the age and metallicity predictions of SPS models 
and their associated errors. 

2. Model ingredients and construction 

Our work focusses on the use of diagnostic 
line indices to determine ages and metallicities for 
SSPs - in particular we will utilize the commonly 
adopted technique of comparing plots of pairs of 
line indices with model grids in our investigation, 
although our results also have important implica- 
tions for other methods, e.g. multiple-index fitting 
or whole-spectrum fitting. 
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Throughout this work we will be using the high 
resolution synthetic s pectra from the BaS TI SPS 
modelH, described in lPercival et al. ( 2009^ . Since 
we have control over all the ingredients of these 
models, we can separately investigate the effects 
of each input element on the resulting spectral en- 
ergy distributions (SEDs) and line indices. 

The two principal ingredients of any SPS model 
are the underlying stellar models, in the form 
of isochrones, which are 'populated' according to 
some initial mass function (IMF) to create each 
stellar population, and the spectral library which 
is used to assign a spectrum to each data point in 
that population. The summing of these individ- 
ual spectra, with appropriate weighting, results in 
the final integrated spectrum (SED) on which di- 
agnostic line strengths, such as H/3 and various 
Fe and other metal lines, can be measured di- 
rectly as equivalent widths (EWs). This was the 
method employed to produce the BaSTI database 
of S EDs, and is used h ere to test the models 
(see IPercival et al.l l2009l for more deta ils) . Al- 
ternat ively, using methods pioneered bv IWorthev 
(|l994l ). some groups use fitti ng functions to deter- 



isochrones are typ ically defined by a f ew hundred 
EPs (e.g. those of lGirardi et"ani200d l. 



mine line strengths (see e.g. iThomas et al.l 12003 ; 



Schiavonll2007t iTrager et all 120081 ). Our study is 



also relevant to these methods since they also rely 
on calibrations from spectral libraries and are con- 
structed in a similar way, i.e. by assigning the rel- 
evant quantity (EW or similar) to each point in a 
population and then summing along the isochrone. 

The choice of isochrones used by SPS modellers 
is undoubtedly a contributory factor to the differ- 
ing predictions from various groups. Underlying 
stellar models (and hence isochrones) from dif- 
ferent sources can vary because of the specifics 
of assumed input physics and differences in stel- 
lar evolution codes. An in-depth discussion of 
these differences is beyond the scope of this pa- 
per, however a detailed comparison of some of the 
m ost commonly used isochrone sets can be found 
in Pietrinferni et al] |2004). One practical dif- 



ference between isochrones from different groups 
is the number of evolutionary points (EPs) that 
define each isochrone - effectively the sampling 
along the isochrone. BaSTI isochrones consist of 
2250 discrete EPs covering all evolutionary stages 
whilst, for example, the commonly used Padova 



Evolutionary points along isochrones are de- 
fined in terms of effective temperature Tgf f , lumi- 
nosity L, and mass Af , from which surface grav- 
ity, \ogg, is also derived. Teff and \ogg, along 
with [Fe/H] (and/or total Z and degree of a- 
enhancement) are used to match, or create by in- 
terpolation, an appropriate spectrum (or fitting 
function) to each EP, since spectra in libraries 
are usually parametrized by these 3 quantities. 
Teff and [Fe/H] (and by implication log;;) must 
be defined on some scale, and there can poten- 
tially be mismatches between the scales arising 
from the underlying stellar models/isochrones and 
those adopted by the spectral library used. This 
is especially relevant when empirical spectral li- 
braries are employed as this requires an evaluation 
of the physical parameters of the observed stars, 
however it is also an issue when using synthetic 
spectral libraries if the stellar models and the at- 
mosphere models are on different scales. Any mis- 
matches in these scales will equally affect SPS 
models using either full SEDs or fitting functions. 

3. Preliminary tests 

As a preliminary evaluation of the typical sys- 
tematic offsets to be expected in effective temper- 
ature scales, we compared the listed Te// values 
for EPs along BaSTI isochrones with those that 
would be 'predi cted' from the empiric al Tef /-color 
calibrations of Alonso et al. I (Il996bt) ffo r dwarfs 
and subdwa rfs) and Alonso et al.l (Il999br) (for gi- 
ants). The lAlonso et all (liggebLligggbh cahbra- 
tions give equations to calculate Te// from vari- 
ous colors, in combination with [Fe/H] - for this 
test the {V — K) color calibration was used, since 
this is the one upon which the atmospheric pa- 
rameters for the MILES empiric al spectral library 



( Sanchez-Blazquez et al. 2006bl used in a later 
Taking the { V — K) color for 



test) are based 

each EP along an isochrone (see lPietrinferni et al 



^http: / /albionc. oa-toramo.inaf.it/ 



|2004 for details of isochrone colors), along with 
[Fe/H], TeffS were calculated from the Alonso 
et al. calibrations and then compared with the 
listed isochrone Te/ /s. Systematic differences were 
found, which have some metallicity dependence, 
and are typically around 70K at solar metallicity 
and up to 160K at [Fe/H] = —2, in the sense that 
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the listed isochrone T^ffS are hotter. These differ- 
ences are similar to the typical quoted systematic 
uncertainties on Teff, which are o f the order 50- 
lOOK (see e.g. lAlonso et al.lll996bl) . 

Ramirez fc Melendej ( 2005 ) determined tem- 



peratures of 135 dwarf and 36 giant FGK stars 
using a method very similar to the Alonso et 
al. studies and found that, in general, the 
agreement in temperature scales was very good, 
and the mean uncertainty in derived temper- 
atures is 75K for dwarfs a nd 60K for giants. 
However, Casagrande et al.l ( |20Q 6) found that 
for 18 stars in common betwee n their study 
and that of Ramirez fc Melendej ( 2005 ). there 
was an average diff e rence of 105±72K, with 
the Casagrande et al. (l2006l ) scale being hotter. 
Casagrande et al.l ( 2006h go on to say 'though 



not negligible, such differences are within the er- 
ror bars of curr ent temperature determinations'. 
More recently Gonzalez Hernandez fc Bonifaciol 
(l2009l) rederived Teff for all stars in the Al onso 
et al. samples ([Alonso et al.l Il996al Il999al) us- 
ing the infrared flux method, and found dif- 
ferences in temperature scales of ~64K and 
~54K for dwar fs and giants respectively, in the 
sense that the Gonzalez Hernandez fc Bonifaciol 
(l2009l) s cale is hotter. However, comparing 



with the iRamirez fc Melendea (|2005[ ) calibration 
they fou nd that for low metalli city dwarfs, the 
[Gonzalez Hernandez fc Bonifacio (2009. ) scale is 
cooler by 

[Fe/H] and logg calibrations are, of course, 
intrinsically linked to Te// when determining 
atmospheric parameters for empirical spectra 
and so systematic uncertainties in T^ff imply 
similar uncertainties on these parameters also. 
Quoted uncertainties are typically of the order 
0.25 dex in logo and 0.15 dex in [Fe/H] (see e.g. 
Soubiran etldll 19981 ). 



In the early stages of creating integrated 
spectra for the BaSTI SPS database we per- 
formed a preliminary comparison between syn- 
thetic spectra from the iMunari et al. Iliooi) spec- 
tral library (as used for the high resolution 
BaSTI SEDs)an d the MILES empiric al spectral 
library jS anchcz- Blazguez et al. 2006b . as used by 
Cenarro et al.l 120081) . For this exercise, we took a 



are from Soubiran et al. ( 19981) (these stars be- 
ing the 'gold standard' for the MILES library). 
To avoid the ambiguity of whether or not to as- 
sume scaled-solar or a-enhanced synthetic spec- 
tra, we used only stars with either [Fe/H]> —0.3 
(assumed to be scaled solar) or [Fe/H]< —1.0 (as- 
sumed to be a-enhanced - note that in the Mu- 
nari library, the level of a-enhancement is fixed at 
[a/Fe]=0.4). By interpolatin g in [Fe/H], T^tf and 
logg amongst spectra in the iMunari et al. [ (120051) 
library, a matching synthetic spectrum was cre- 
ated for each star in the MILES subsample. Line 
strengths were then measured directly on the two 
sets of spectra using t he LECTOR program of A. 
Vazdekislfl (as done in Percival et al. 20091 ) and a 
mean offset found of ~ 0.4 dex in H/3, the prin- 
cipal age indicator for SSPs, in the sense that 
the MILES spectra have larger values. Significant 
mean offsets were also found in the main metallic- 
ity indicators, e.g. AMgb ~ 0.6, A[MgFe] ~ 0.5 
and AFe5406 ~ 0.2, in the sense that the MILES 
spectra have smaller values. One possible interpre- 
tation is that these systematic offsets are related to 
the calibration of atmospheric parameters, which 
may differ between the two spectral libraries. 
The systematic offsets that w e find are qualita- 
tively similar to those found bv lMartins fc Coelhd 
(|2007l ). who performed a comparison of spectral 
indices measured on several theoretical and em- 
pirical spectral lib raries, using a slightly dif ferent 



method to ours. iMartins fc Coelhd (|2007l) gen 



subset of the MILES spectra with zero redden- 
ing {E{B ~ V) — 0.0) and with the designation 
SKC, indicating that the atmospheric parameters 



erally found that disagreements between spectral 
libraries have some temperature dependence, and 
that the largest offsets often occur for the cooler 
stars (defined as Tg// < 4500K in their work). 
Exploring the impact that these differences can 
have on SPS models is a key motivation for the 
work presented in this paper. 

In the following tests we examine the effects 
of altering each of the atmospheric parameters in 
the the BaSTI SPS models separately, within the 
typical uncertainties, whilst holding all other ele- 
ments constant - these uncertainties are taken to 
be iblOOK in Te//, ±0.25 dex in \ogg and ±0.15 
dex in [Fe/H]. We stress that these tests are ex- 
ploring the effects of a mismatch in scales between 
the underlying stellar models (and isochrones) and 
the adopted spectral hbrary in the SPS models. 



^see [http: / /www.iac.es/galeria/vazdekis/| 
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which potentially give rise to systematic errors in 
derived ages and metallicities for SSPs. When 
presenting our results we are making the implicit 
assumption that the differential behavior we find 
would be quantitatively the same, or very similar, 
for other SPS models, whether based on full inte- 
grated spectra or on fitting functions methods to 
determine line strengths. 

4. The tests 

All the tests described below were performed 
on two test-case SSPs with ages t=lA Gyr, to rep- 
resent a typical old elliptical galaxy, and t=A Gyr, 
which is representative of the intermediate ages 
found in th e sam ple of ellipticals studied by 
Trager et al. I 12000) . In both cases the scaled- 
solar abundance, solar metallicity models were 
used {Z = 0.0198 and [Fe/H] = -f0.06) with no 
convective overshooting, and with the Reimers 
mass-loss parameter, rj, set at 0.2. High resolu- 
tion (lA/pixel) integrated spectra for the refer- 
ence models were created, as described in Percival 
et al. (2009), where the interested reader can also 
find more details on the underlying stellar models. 
We note here that we have also applied all the fol- 
lowing tests to low metallicity a-enhanced SSPs, 
typical of Galactic globular clusters, and obtain 
quantitatively very similar results. 

When constructing the integrated spectrum for 
each test, the underlying isochrone is left un- 
changed, but individual spectra are assigned to 
each EP along the isochrone with the appropriate 
offset in each atmospheric parameter, as listed be- 
low. Each atmospheric parameter was tested for 
two instances - one in which the whole isochrone 
is affected (i.e. stars at all evolutionary phases), 
and one in which only the giant stars are affected. 
For these purposes we took all EPs with, simul- 
taneously, Teff <5000K and logg <3 to represent 
giants. The justification for treating giants only 
as a separate case is that, for empirical T^/f cali- 
brations, dwarfs and giants are often treated sepa- 
rately and consequently have separate Te//-color 
relationships. It is also generally the case that 
cool stars are harder to model theoretically due to 
the effects of phenomena such as molecular opac- 
ities, convection and mass loss which become in- 
creasingly important for cooler stars, hence there 
is more likely to be a discrepancy in temperature 



scales for cool RGB and AGB stars. Tests were 
performed on the two reference SSPs as follows: 

Test 1 - Spectra assigned with Tf,ff increased 
by lOOK compared to isochrone T^ffS (whole 
isochrone). 

Test 2 - Spectra assigned with T^ff increased 
by lOOK compared to isochrone TeffS, for giants 
only (i.e. Tejf <5000K and logg <3 on isochrone). 

Note that we could not perform fully equivalent 
tests in which Te// was decreased by lOOK. This 
is because there is a low temperature cut off in 
the spectral library use d fo r the BaSTI SPS m od- 



els ( Munari et al. 2005 , see Percival et al. 20091 for 



more details) which means that the coolest stars 
could not be consistently included in the test as 
an offset of -lOOK would take them below this 
threshold. However, we did perform this test us- 
ing a truncated version of the 4 Gyr SSP, i.e. with 
a cut at lOOK above the threshold T^/f, so that 
this test could be done by assigning spectra with 
both increased and decreased Tefj values. 

Tests 3a, b - logg increased/decreased by 0.25 
dex (whole isochrone) 

Tests 4a, b - logg increased/decreased by 0.25 
dex (giants only) 

Tests 5a, b - [Fe/H] increased/decreased by 0.15 
dex (whole isochrone) 

Tests 6a, b - [Fe/H] increased/ decreased by 0.15 
dex (giants only) 

Two further tests were also performed in which 
two 'worst-case' scenarios were created. For the 
first of these, two parameters were altered simul- 
taneously over the whole isochrone, namely T^fj 
increased by lOOK and logg decreased by 0.25 dex. 
For the second of these tests all three parameters 
were altered - Tef / was increased, logg decreased 
and [Fe/H] increased by 0.15 dex. The purpose of 
these tests was to check whether the effects on line 
indices are additive when two or three atmospheric 
parameters are offset simultaneously. 

Line strengths for 23 diagnostic indices were 
then measured for the reference cases and test 
cases in the same way as before, i.e. by measuring 
equivalent widths directly on the integrated spec- 
tra themselves. Indices use d are the 21 i ndices 
defined by the bandpasses in Trager et al. ( 19981 ) 



(noting that, by convention, CNi, CN2, Mgi, Mg2, 
TiOi and Ti02 are quoted in magnitudes rather 
than as EWs in Angstroms) plus HSp and H'^f as 
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defined in lWorthev fc Ottavianil (jlQQTT) . We stress 
that all numbers quoted are line strengths as mea- 
sured directly on the integrated spectra and are 
not transformed onto the Lick system. 

5. Results 

The results of the tests on the 4 Gyr and 14 Gyr 
populations are tabulated in Tables [1] and [2] re- 
spectively, where the second column of each ta- 
ble shows the measured index strengths for the 
unaltered reference SSP, and subsequent columns 
list the offset in each index between the reference 
SSP and each of the spectra resulting from tests 
1-6. The quantities listed are offsets which must 
be added to the reference values to reproduce the 
results of each test. 

5.1. Effects of altering T^ff 

Results from tests 1 and 2, i.e. T^ff increased 
for the whole isochrone and for giants only, respec- 
tively, are shown in the third and fourth columns 
of Tables [1] and [21 For both test-case SSPs it can 
be seen that H/3 increases significantly when T^f / 
is increased for the whole isochrone, as would be 
expected. However when Te// is increased for gi- 
ants only, the H/3 value increases slightly for the 
14 Gyr SSP and actually decreases for the 4 Gyr 
SSP - we attribute this to the fact that for giant 
stars at these low temperatures {T^ff < 5000K) 
and gravities {logg <3), increasing T^fj by lOOK 
increases the continuum level around the H/3 line 
but does not significantly increase the strength of 
the line itself, and so the measured H/3 value can 
decrease. The higher order Balmer lines HSp and 
HjE behave in a qualitatively similar way to H/3 
in that they increase significantly when T^f / is in- 
creased for the whole isochrone, but increase only 
marginally, or decrease slightly, for the giants-only 
case. 

For all but one of the metal indicators (in- 
cluding Mg5 and all the Fe lines) index strengths 
decrease when T^fj is increased for the whole 
isochrone, the exception being C24668 for the 
4 Gyr SSP, which increases slightly. When T^ff 
is increased for giants only, offsets are generally 
smaller, with the exception of C24668, but for 
many indices offsets go in the opposite sense to 
those for the whole-isochrone case, e.g. CNl, CN2 
and most of the Fe lines. 



Results of the decreased T^fj test (using the 
truncated 4 Gyr SSP) are not listed, as the mea- 
sured values of offsets would not be applicable to 
the full SSP. However we note that, for all 23 mea- 
sured indices, offsets for the decreased Te// trun- 
cated SSP were found to be symmetrically oppo- 
site with respect to the increased Tg// truncated 
SSP case, implying that this would also be the 
case for the full SSP. 

To demonstrate the effect of these offsets on es- 
timated ages and metallicities of SSPs, we have 
constructed representative index-index grids from 
the BaSTl models for a range of ages and metallici- 
ties and overplotted the test-case SSPs - FiguresIU 
m and E] show grids of H/3 vs. Mg6, Fe5406 and 
C24668, respectively. In each figure the left hand 
panel shows the results of the whole-isochrone 
tests and the right hand panel shows the results of 
the giants-only tests. It can be seen that increas- 
ing Teff by lOOK for the whole isochrone makes 
the apparent age of the 4 Gyr SSP approximately 
0.75 Gyr younger whilst the 14 Gyr SSP looks 
2.5—3 Gyr younger, an apparent reduction in age 
of ~20% for both populations. The effect on ap- 
parent ages for the giants-only test is small, alter- 
ing absolute ages by no more than ~3%, however 
since the offsets for the 4 Gyr and 14 Gyr popula- 
tions go in opposite senses, the relative ages also 
change. 

Assessing the effect of temperature changes on 
derived metallicities is more complicated, since 
it depends on which index is used and whether 
metallicity is defined in terms of [Fe/H] or total 
metallicity, Z. The three grids displayed in Fig- 
ures [U [2] and [3] demonstrate the problem. 

For a lOOK increase in T^ff, the inferred [Fe/H] 
from the H/3-Mg6 plane is ~ 0.1 dex lower than 
the reference SSP, from the H/3-Fe5406 plane it 
is virtually unchanged, and from the H/3-C24668 
plane [Fe/H] apparently increases by ^ 0.05 dex. 
Offsets are similar for both the 'whole-isochrone' 
and 'giants-only' cases. The H/3-C24668 grid also 
demonstrates another important point - for the 
14 Gyr SSP, the offset for the increased lOOK case 
decreases the absolute value of C24668 slightly, 
however the inferred [Fe/H] increases as this grid 
is not orthogonal (as is also the case for all other 
commonly used age /metallicity diagnostic grids). 
Similarly constructed grids (not displayed here) 
show that, with increasing T^ff, Ca4227 behaves 
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Fig. 1.— H/3 vs. Mg6 grid for BaSTI scaled-solar 
SSP models for ages 3, 6, 8, 10 and 14 Gyr (age 
increasing from the top downwards) and [Fe/H] = 
-0.66, -0.35, +0.06 (solar) and +0.40 (increas- 
ing from left to right). The left-hand panel shows 
results of the whole-isochrone tests and the right- 
hand panel shows results of the giants-only tests. 
The two test reference SSPs (solar metallicity, 
4 Gyr and 14 Gyr) are marked as solid squares. 
Results of various tests are labelled on the dia- 
grams (see text for details). 
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Fig. 2. — H/3 vs. Fe5406. Panels, grid points and 
symbols are the same as for Figure [T] 
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Fig. 3. — H/3 vs. C24668. Panels, grid points and 
symbols are the same as for Figured] 
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[MgFe] 

Fig. 4. — Fe5406 vs. [MgFe] diagram for 
14 Gyr SSP scaled solar (solid line) and a- 
enhanced (dashed line) models, joined at approx- 
imately equal [Fe/H] (dotted lines) for [Fe/H] ~ 
+0.06, —0.3, —0.7, —1.0 (decreasing from the top 
downwards). Test reference SSPs are marked as 
squares, and other symbols are as for Figure [T] 
for the 4 Gyr (open symbols) and 14 Gyr (solid 
symbols) cases. 
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in the same way as Mg6, [MgFe] behaves as Fe5406 
(but traces total Z rather than [Fe/H] - see 
IPercival et aD l2009f l and CNl and CN2 behave 
similarly to C24668. The issue of the contradic- 
tory behavior of various diagnostic metal lines and 
the impact on inferred chemical composition will 
be discussed further in Section [71 

5.2. Effects of altering logg 

It was found that, for the logy tests, offsets in 
all indices were found to be symmetrical for the 
increased/decreased cases with respect to the ref- 
erence case, hence the results are listed in Tables [T] 
and m as single values for each pairing of tests 
3a/3b and 4a/4b, with the appropriate ± or =F 
sign. For both test-case SSPs it can be seen that 
H/3, and H7i^^ increase significantly when logg 
is decreased by 0.25 dex (and decrease when logg 
is increased), affecting derived ages in a quantita- 
tively similar way to a lOOK increase in Te/ / . Mg6, 
Ca4227, Fe5406, [MgFe] and C24668 also behave in 
a quantitatively very similar way with decreasing 
logg to the increased T^j f case. Figures [U [2] and 
[3] show that, with increasing logg, [Fe/H] inferred 
from the H/?~Mg6 plane increases by ^ 0.1 dex, 
whilst that from the H/3-[Fe5406] plane is virtually 
unchanged and from the H/3-C24668 plane [Fe/H] 
apparently decreases by ~ 0.1 dex. 

5.3. Effects of altering [Fe/H] 

As for the \ogg and Te// tests, the [Fe/H] tests 
produced offsets in all indices that were found to 
be symmetrical for the increased/decreased cases 
with respect to the reference case. Results are 
listed in Tables [T] and [2] as single values for each 
pairing of tests 5a/5b and 6a/6b, with the ap- 
propriate ± or =F sign. As would be expected, 
all metal lines increase in strength when [Fe/H] is 
increased and decrease when [Fe/H] is decreased. 
For the whole-isochrone tests, the inferred change 
in [Fe/H] from the diagnostic diagrams is con- 
sistent with the modelled change in [Fe/H], i.e. 
~ 0.15 dex, whilst for the giants-only test the in- 
ferred change in [Fc/H] is generally about half this 
level. 

An important point to notice is that H/3 in- 
creases significantly when [Fe/H] is increased over 
the whole isochrone (and decreases when [Fe/H] 
is decreased) inducing ~ 12% change in derived 



ages for a 0.15 dex change in [Fe/H], however for 
the giants-only case H/3 is virtually unchanged. 
Rather puzzlingly, HSp and Hjf behave in the 
opposite sense to H/3 in that they both decrease 
significantly with increased [Fe/H], and vice versa. 
Hdp and H'Jf also display quantitatively simi- 
lar offsets for the whole-isochrone and giants-only 
cases, unlike the situation for H/3. The implica- 
tions of this contradictory behavior will be dis- 
cussed further in Section [71 

5.4. 'Worst-case' scenarios 

The results from the 'worst-case tests', in which 
two or three parameters were altered simultane- 
ously, are not listed - however in both cases the 
offsets in all indices were found to be additive. For 
the two parameter test (altering T^f / and \ogg si- 
multaneously) resultant offsets are simply the sum 
of the offsets found for tests 1 and 3b, whilst for 
the three parameter test {T^ff, logg and [Fe/H]) 
offsets are the sum of those from tests 1, 3b and 
5a. Results of the two 'worst-case' tests are also 
illustrated in Figures [H [2] and [3l 

6. Significance of results: comparison with 
typical observational errors 

In order to assess the significance of our results 
for observed stellar populations, we briefly com- 
pare the line index offsets detailed above with typ- 
ical observational errors for extragalactic systems. 
Table [3] shows the mean observational errors on 
various line index strengths fr om three sources. 
Sanchez-Blazquez et al. (|2006a[ ) present a study 



of 98 early type (E and SO) galaxies in the lo- 
cal field, small groups, and some Virgo and Coma 
cluster members, whilst [Trager et aL [ (120081 ) study 
12 elliptical and SO galaxies in the Coma clus- 
ter. For hig her redshift s ysteni s we used data from 
[Sanchez-Blazquez et al. ( 20091 ) . which comprises a 
catalogue of 215 red sequence galaxies in clusters 
and groups, with rcdshifts between z ~ 0.45 and 
z ~ 0.75. For the ISanchez-Blazguez et al" ( 2009[ ) 
sample, observational errors on all quoted indices 
are slightly larger than the largest offsets found 
in our work, so that the systematic effects that 
we find would be a significant, but not the domi- 
nant, source of error for these systems. However, 
for the more local systems observational errors for 
all the key diagnostic indices are significantly less 
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than the offsets caused by the mismatches in stel- 
lar parameters detailed above. Most importantly, 
the observational errors on HSp, Hjf and H/3 are 
smaller by a factor of ~ 2 than the offset induced 
by a lOOK mismatch on temperature scales, which 
would lead to a systematic shift in the inferred 
ages for these systems. Similarly, observational 
errors on Ca4227, Mg6 and most of the Fe lines 
are significantly smaller than the offsets in indices 
caused by a mismatch in logg or [Fe/H] scales. It 
should be borne in mind that any offsets resulting 
from stellar parameter mismatches induce system- 
atic rather than random errors. 

We remind the reader that all our results dis- 
cussed so far pertain to intermediate age and old 
stellar populations. At the suggestion of the ref- 
eree we also performed several tests on a 500 Myr, 
solar metallicity population. At this age (and 
younger) stellar populations do not have a red gi- 
ant branch as such, and so we performed the basic 
tests of increasing Teff by lOOK, decreasing \ogg 
by 0.25 dex and altering [Fc/H] by 0.15 dex, for 
the whole isochrone case only. The offsets in in- 
dices found for this younger age population are 
generally a factor ^ 2 smaller than those found 
for the 4 Gyr and 14 Gyr SSPs and so in most 
cases are less significant than the observational er- 
rors quoted in Tabled Also, the strength of the 
Balmer lines increases much more rapidly with de- 
creasing age for ages below ~ 1 Gyr and so even 
a substantial change in H/3 of 0.2 dex (larger than 
any of the offsets found here) only changes the in- 
ferred age by 10 — 20 Myr. 

6.1. Effect on broad-band colors 

The aim of this paper is primarily to investi- 
gate the effect of mismatches in stellar parameters 
on diagnostic line indices, however we have also 
made a brief study of the effect on broad-band 
colors for the main tests outlined in Section [D In 
general, a mismatch in the \ogg scale has a negligi- 
ble effect on colors whilst the effects of a 0.15 dex 
mismatch in [Fe/H] are small and are typically 
within the likely photometric errors, i.e. < 0.02 
mag in all the optical and near-infrared bands. As 
might be expected, the effect of increasing T^ff 
by lOOK (Test 1 above) is to make all the broad- 
band colors bluer. For the intermediate age and 
old populations studied here, the effect is small 
in the the optical bands, ~ 0.01 mag in {U — B) 



and {B — V), and 0.02 mag in {R — I), imply- 
ing a minimal change in the inferred age, however 
the effect becomes greater at longer wavelengths, 
making the (J — K) color bluer by ~ 0.04 mag. 
Hence broader baseline colors such as {V — K) 
and {B — K) are affected more strongly by an in- 
crease in Teff, the colors becoming bluer by at 
least 0.1 mag. As a result, a lOOK increase in Teff 
can significantly decrease the inferred metallicity, 
rather than making a population look younger, if 
a color-color diagram such as {V — I)/{V — K) 
or {B — K)/(J — K) is used (see Fi gure 1 of. re- 
spect ively, ISalaris &: CassisilboOTi and James et al.l 
20061 ). 

7. Summary and discussion 

We have investigated a potentially significant 
source of systematic error which can affect the 
ages and metallicities of stellar populations de- 
rived from SPS modelling, caused by systematic 
uncertainties in the three principal stellar atmo- 
spheric parameters Teff, logg and [Fe/H]. In prac- 
tice, we have tested the effects of a mismatch 
in temperature, metallicity and gravity scales be- 
tween the underlying isochrones and the spectra 
(or fitting functions) used to construct SPS mod- 
els. We have done this by constructing integrated 
spectra for two solar metallicity SSP models aged 
4 Gyr and 14 Gyr, and applying offsets of lOOK 
in Teff, 0.25 dex in logg and 0.15 dex in [Fe/H], 
which we take as typical zero-point uncertainties 
in these key atmospheric parameters. We have 
then quantified the effect on various diagnostic line 
indices by measuring EWs on the resultant spec- 
tra and calculating offsets between the unaltered 
reference SSPs and the altered versions. 

We note here that the magnitude of the offsets 
we find as a result of altering the stellar parame- 
ters for SSPs, as described, is quantitatively very 
similar to the differences in line index strengths 
found between the various spectral libraries, de- 
tailed in Section [3] Hence our results are con- 
sistent with the differences between the spectral 
libraries being largely caused by zero-point differ- 
ences in their Teff, logg or [Fe/H] scales. We take 
this as an indication that our results give a good 
estimate of the likely systematic errors in line in- 
dex strengths inherent in any population synthesis 
model, due to possible zero-point mismatches in 
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stellar atmospheric parameter scales. 

Using simple index-index diagrams to make a 
preliminary assessment of the impact on stehar 
population parameters inferred from the models, 
we find that absolute ages derived from the H/3 in- 
dex can easily be affected at the 20% level for both 
old and intermediate age populations. Relative 
ages can also be affected, albeit at a lower level. 
Inferred ages are more complicated to interpret if 
the H6f and H7i? indices are also considered be- 
cause of their behavior in response to a systematic 
shift in [Fe/H], which goes in the opposite sense 
to that of H/3. The inferred systematic errors in 
[Fe/H] and/or Z are also hard to quantify because 
of the opposite behavior of certain key metal in- 
dicators in response to systematic offsets in Te// 
and \ogg. 

This behavior has implications for methods 
which fit simultaneously to several indices (or per- 
form a full SED fit) to derive ages and metallicities 
of stellar populations, since a failure to fit several 
indices simultaneously could, spuriously, be inter- 
preted as an indication of non-solar abundance ra- 
tios. Also, mismatches between observational data 
and model SSPs are often taken as evidence for 
the presence of components which are not fully 
accounted for in the models, such as missing (or 
extreme) stellar evolutionary stages or a compos- 
ite population. 

Percival et al. (I2OO9I) demonstrated that, for 
SSPs, the Fe5406 index traces Fe only, whilst 
[M gFe] traces total met allicity, Z (as first noted 



bv iThomas et al. 2003[ ). These two indices, in 
combination therefore provide an estimate of the 
level of a-enhancement, or simply whether a pop- 
ulation has non-solar abundance ratios. Figured 
shows the Fe5406-[MgFe] plane with lines from 
the BaSTl 14 Gyr (constant age) scaled-solar and 
a-enhanced models, joined at points of approx- 
imately equal [Fe/H] (note that at constant Z, 
[Fe/H] is lower for a-enhanced models than the 
corresponding scaled-solar ones). Overplotted are 
the results from the whole-isochrone tests for in- 
creased Teff (lOGK), increased \ogg (0.25 dex) and 
increased/decreased [Fe/H] (0.15 dex). In this dia- 
gram one expects that any deviation in abundance 
ratios will move points horizontally, i.e. any degree 
of a-enhancement moves points from the scaled- 
solar line on the left, towards the g-enhanced line 
on the right (see Figure 9 of lPercival eral]l2009l 



- note that lines of different age are completely 
degenerate in this diagram). Here it can be seen 
that altering any of the atmospheric parameters 
simply moves the scaled-solar SSP points along 
the scaled-solar line, changing the inferred [Fe/H] 
but not altering the inferred abundance ratios for 
other elements to non-solar ratios. This is an im- 
portant point to notice since individual element 
abundances, including several a elements, appear 
to alter substantially, as demonstrated by the mea- 
sured offsets in the various line indices listed in 
Tables [I] and H 

In conclusion, we urge caution against the 
over-interpretation of stellar population param- 
eters from line index data, in terms of the in- 
ferred scaled-solar or non-scaled-solar abundance 
ratios and also the inferred presence of a compos- 
ite population, especially when multiple-index or 
fuU-SED fitting methods are employed. We find 
that, for SSPs, Fe5406 in combination with [MgFe] 
provides the most robust indication of non-solar 
abundance ratios. We remind the reader that our 
results potentially impact on all SPS methods, 
whether fitting functions or full SEDs are em- 
ployed. Measured offsets for 23 commonly used 
diagnostic line indices are provided, and we en- 
courage the user to determine the overall impact 
on their observational data and preferred fitting 
method. 

We thank the anonymous referee for a con- 
structive report and some useful suggestions which 
helped to put our results in context. S.M.P. 
would like to express heartfelt thanks to Elaine 
Smith-Freeman for many useful discussions and 
for providing the initial motivation to do this work. 
S.M.P. acknowledges financial support from the 
Science & Technology Facilities Council (STFC) 
through a Postdoctoral Research Fellowship. 
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Table 1 

Results of tests on the 4Gyr SSP. The second column lists EWs measured on the 

reference ssp spectrum and subsequent columns list the offset in each index between 
the reference spectrum and each of the test spectra (see text for details). 



Te//+100K logc;±0.25 [Fe/H]±0.15 

Full iso Giants Full iso Giants Full iso Giants 
Index Ref (1) (2) (3a,b) (4a,b) (5a,b) (6a,b) 





0. 


.083 


0. 


.228 


-0. 


.109 


TO. 


.098 


TO, 


.011 


TO. 


.067 


TO, 


.081 


H7f 


-1, 


.208 


0, 


.315 


-0, 


.085 


TO, 


.181 


TO, 


.031 


TO, 


.078 


TO, 


.079 


CNi 


0, 


.002 


-0, 


.003 


0. 


.007 


TO. 


.001 


TO, 


.001 


±0, 


.007 


±0, 


.005 


CN2 


0, 


.058 


-0, 


.004 


0. 


.005 


±0. 


.000 


±0, 


.000 


±0. 


.009 


±0, 


.005 


Ca4227 


1, 


.049 


-0, 


.141 


-0, 


.098 


±0. 


.134 


±0, 


.069 


±0. 


.135 


±0, 


.109 


G4300 


7, 


.198 


-0, 


.181 


0, 


.126 


±0. 


.036 


TO, 


.033 


±0, 


.106 


±0, 


.004 


Fe4383 


5, 


.743 


-0, 


.272 


0. 


.029 


±0. 


.185 


±0, 


.048 


±0. 


.243 


±0, 


.094 


Ca4455 


1, 


.747 


-0, 


.063 


0, 


.004 


±0. 


.014 


±0, 


.010 


±0, 


.128 


±0, 


.062 


Fc4531 


4, 


.696 


-0, 


.064 


0, 


.054 


TO. 


.043 


TO, 


.036 


±0, 


.308 


±0, 


.145 


C24668 


2, 


.632 


0. 


.008 


0. 


.170 


TO. 


.187 


TO, 


.168 


±0. 


.512 


±0, 


.311 


H/3 


2, 


.015 


0, 


.130 


-0, 


.041 


TO, 


.078 


TO, 


.025 


±0, 


.087 


±0, 


.009 


Fe5015 


7, 


.005 


-0, 


.045 


0, 


.108 


TO, 


.141 


TO, 


.115 


±0, 


.558 


±0, 


.272 


Mgi 


0, 


.103 


-0, 


.006 


-0. 


.001 


±0, 


.003 


±0, 


.002 


±0. 


.011 


±0, 


.008 


Mg2 


0, 


.252 


-0, 


.019 


-0, 


.010 


±0, 


.019 


±0, 


.009 


±0. 


.021 


±0, 


.012 


Mg6 


3, 


.755 


-0, 


.323 


-0. 


.194 


±0, 


.482 


±0, 


.206 


±0, 


.277 


±0, 


.145 


Fe5270 


3, 


.524 


-0, 


.074 


0. 


.015 


±0. 


.041 


TO, 


.005 


±0. 


.272 


±0, 


.122 


Fe5335 


4, 


.081 


-0, 


.081 


0, 


.015 


±0, 


.045 


TO, 


.006 


±0, 


.351 


±0, 


.192 


Fc5406 


2, 


.188 


-0, 


.086 


-0, 


.022 


±0, 


.052 


±0, 


.015 


±0, 


.169 


±0, 


.103 


Fe5709 


1, 


.388 


-0, 


.014 


0. 


.026 


TO, 


.031 


TO, 


.027 


±0. 


.134 


±0, 


.075 


Fe5782 


0, 


.879 


-0, 


.012 


0, 


.016 


TO, 


.013 


TO, 


.013 


±0, 


.098 


±0, 


.061 


Na D 


2, 


.813 


-0, 


.194 


-0, 


.106 


±0, 


.178 


±0, 


.071 


±0, 


.265 


±0, 


.199 


TiOi 


0, 


.019 


-0, 


.005 


-0. 


.004 


±0, 


.001 


±0, 


.001 


±0. 


.002 


±0, 


.002 


Ti02 


0, 


.049 


-0, 


.007 


-0. 


.006 


±0, 


.000 


±0, 


.000 


±0. 


.005 


±0, 


.004 
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Table 2 

Results of tests on the 14Gyr SSP. Columns are as for Table [U 



Te//+100K logg±0.25 [Fe/H]±0.15 

Full iso Giants Full iso Giants Full iso Giants 
Index Ref (1) (2) (3a,b) (4a,b) (5a,b) (6a,b) 





-1, 


.755 


0. 


.282 


-0. 


.026 


TO. 


.083 


±0, 


,003 


TO, 


.152 


TO, 


,102 




-3 


,151 


0. 


.273 


0. 


.021 


TO. 


.132 


TO, 


,018 


TO, 


.101 


TO, 


,093 


CNi 


0, 


.056 


-0, 


.007 


0. 


.005 


TO. 


.004 


TO. 


,001 


±0, 


.013 


±0, 


,005 


CN2 


0, 


.114 


-0, 


.010 


0. 


.003 


TO. 


.002 


TO. 


,001 


±0. 


.017 


±0, 


,006 


Ca4227 


1, 


.672 


-0, 


,271 


-0. 


.144 


±0. 


.234 


±0, 


,073 


±0. 


.200 


±0, 


,134 


G4300 


8, 


.287 


-0, 


.035 


0. 


.069 


TO. 


.040 


TO, 


,042 


±0. 


.015 


±0, 


,008 


Fe4383 


7, 


.553 


-0, 


,338 


-0. 


.061 


±0. 


.198 


±0, 


,029 


±0. 


.263 


±0, 


,095 


Ca4455 


2, 


.167 


-0, 


,101 


-0. 


.017 


±0. 


.026 


±0, 


,008 


±0. 


.150 


±0, 


,062 


Fe4531 


5, 


.319 


-0, 


,098 


0. 


.025 


TO. 


.033 


TO, 


,041 


±0. 


.322 


±0, 


,143 


C24668 


3, 


.228 


-0, 


,009 


0. 


.143 


TO. 


.258 


TO, 


,168 


±0. 


.582 


±0, 


,301 


H/3 


1, 


.135 


0. 


.134 


0. 


.010 


TO. 


.072 


TO, 


,023 


±0. 


.067 


±0, 


,008 


Fe5015 


7, 


.642 


-0, 


.066 


0. 


.077 


TO. 


.150 


TO, 


,116 


±0. 


.569 


±0, 


,265 


Mgi 


0, 


.144 


-0, 


.012 


-0. 


.004 


±0. 


.006 


±0, 


,002 


±0. 


.014 


±0, 


,007 


Mg2 


0, 


.338 


-0, 


.029 


-0. 


.015 


±0. 


.026 


±0, 


,009 


±0. 


.025 


±0, 


,012 


Mg6 


4, 


.981 


-0, 


.421 


-0. 


.277 


±0. 


.545 


±0, 


,192 


±0. 


.329 


±0, 


,145 


Fe5270 


4, 


.107 


-0, 


,112 


-0. 


.014 


±0. 


.063 


TO, 


,009 


±0. 


.273 


±0, 


,112 


Fe5335 


4, 


.680 


-0, 


,117 


-0. 


.012 


±0. 


.052 


TO, 


,011 


±0. 


.374 


±0, 


,174 


Fe5406 


2, 


.640 


-0, 


,128 


-0. 


.043 


±0. 


.069 


±0, 


,012 


±0. 


.191 


±0, 


,096 


Fe5709 


1, 


.557 


-0, 


,021 


0. 


.020 


TO. 


.039 


TO, 


,026 


±0, 


.139 


±0, 


,068 


Fe5782 


1, 


.019 


-0, 


,021 


0. 


.012 


TO. 


.016 


TO, 


,014 


±0. 


.107 


±0, 


,055 


Na D 


3, 


.711 


-0, 


,306 


-0. 


.160 


±0. 


.269 


±0, 


,067 


±0. 


.317 


±0, 


,190 


TiOi 


0, 


.022 


-0, 


,006 


-0. 


.004 


±0. 


.001 


±0, 


,001 


±0. 


.002 


±0, 


,002 


Ti02 


0, 


.056 


-0, 


,009 


-0. 


.006 


±0. 


.000 


TO, 


,000 


±0. 


.006 


±0, 


,004 
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Table 3 

Mean o bservati onal errors on line-strengths from 3 sour ces: ISanchez-Blazquez et alI 
(|2006aI ') (SB06). IsAnchez-Blazquez et alI (|2009l ) (SB09) and ITrager et al.I (|2008[ ) f T08) 



Index 


SB06 


SB09 


T08 


RSp 


0.107 


0.303 


0.100 


H7f 


0.108 


0.314 


0.111 


CNi 


0.006 




0.004 


CN2 


0.007 


0.016 


0.004 


Ca4227 


0.076 


0.250 


0.067 


G4300 


0.184 


0.438 


0.128 


Fe4383 


0.208 


0.721 


0.152 


Ca4455 


0.093 


0.396 


0.086 


Fe4531 


0.131 


0.613 


0.107 


C24668 


0.253 




0.165 


H/3 


0.079 




0.066 


Fe5015 


0.189 




0.142 


Mg6 


0.143 




0.068 


Fe5270 






0.076 


Fe5335 






0.085 



Bruzual, A. G. 2007, in lAU Symposium, Vol. 
241, lAU Symposimii, ed. A. Vazdekis & R. F. 
Peletier, 125-132 

Bruzual, G. & Chariot, S. 2003, MNRAS, 344, 
1000 

Carter, D., Smith, D. J. B., Percival, S. M., 
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